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Abstract 

Using the recently introduced method to calculate bubble abundances in an eternally inflating 
spacetime, we investigate the volume distribution for the cosmological constant A in the context 
of the Bousso-Polchinski landscape model. We find that the resulting distribution has a staggered 
appearance which is in sharp contrast to the heuristically expected flat distribution. Previous 
successful predictions for the observed value of A have hinged on the assumption of a flat volume 
distribution. To reconcile our staggered distribution with observations for A, the BP model would 
have to produce a huge number of vacua in the anthropic range AA^ of A, so that the distribution 
could conceivably become smooth after averaging over some suitable scale 5 A <C AA^. 
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I. INTRODUCTION 



The cosmological constant problem is one of the most intriguing mysteries that we now 
face in theoretical physics. The observed value of the cosmological constant A is many orders 
of magnitude smaller than theoretical expectations and is surprisingly close to the present 
matter density of the universe, 1 

A ~ Pm0 ~ 1(T 120 . (1) 

As of now, the only plausible explanation for these enigmatic facts has been given in terms of 
the multiverse picture, which assumes that A is a variable parameter taking different values 
in different parts of the universe 0, H, IE HL IM IM IE H( • The probability for a randomly picked 
observer to measure a given value of A can then be expressed as 

P o6s (A) ocP(A)n o6s (A), (2) 

where P(A) has the meaning of the volume fraction of the regions with a given value of A 
and n b s (A) is the number of observers per unit volume. 2 Disregarding the possible variation 
of other "constants" and assuming that the density of observers is roughly proportional to 
the fraction of matter clustered in large galaxies, 

n obs (A) ex f G (A), (3) 

one finds that the function n obs (A) is narrowly peaked around A = 0, with a width 

AA A ~ 100A ~ 10 -118 . (4) 

In general, the volume factor P(A) depends on the unknown details of the fundamental 
theory and on the dynamics of eternal inflation. However, it has been argued jE Ell ^ na ^ ^ 
should be well approximated by a flat distribution, 

P(A) « const. (5) 

The reason is that the anthropic range (@J), where the function n b s (A) is substantially 
different from zero, is much narrower than the full range of variation of A, which is expected 
to be set by the Planck scale. A smooth function varying on this large characteristic scale 
will be nearly constant within the tiny anthropic interval. 
Combination of Eqs. ©-© yields the distribution 

P obs (A) cx f G (A), (6) 

which can be readily calculated using the Press-Schechter approximation for fa. The ob- 
served value of A is within the 2a range of this distribution - an impressive success of the 
multiverse paradigm. One should keep in mind, however, that the successful prediction for 



1 Here and below we use reduced Planck units, Mp/87r = 1, where M p is the Planck mass. 

2 .P(A) is often called the prior probability. Here we avoid this terminology, since it is usually used to 
characterize one's ignorance or prejudice, while the volume factor P(A) should be calculable, at least in 
principle. 
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A hinges on the assumption of a flat volume distribution (J3J). We emphasize that the form of 
the volume distribution is important. If, for example, one uses P(A) oc A instead of (0), the 
2cr prediction would be 10 < A/A < 500 and the observed value of A would be ruled out 
at a 99.9% confidence level [llj]. The heuristic argument for a flat distribution (JHJ) sounds 
plausible, but it needs to be verified in specific models. 

The simplest model with a variable effective cosmological constant is that of a scalar field 
with a very slowly varying potential V{<j>) HQ In such models, A takes a continuum 
range of values. It has been verified that the resulting volume distribution for A is indeed 
flat for a wide range of potentials (l3l Q, ^| . The main challenge one has to face in this 
type of model is to justify the exceedingly flat potential which is required to keep the field 
(f) from rolling downhill on the present Hubble time scale. 

A model with a discrete spectrum of A was first suggested by Abbott |l6| . He considered a 
scalar field with a "washboard" potential, having many local minima separated by barriers. 
Transitions between the minima could occur through bubble nucleation. An alternative 
discrete model, first introduced by Brown and Teitelboim [lTj , assumes that the cosmological 
constant can be expressed as 

A = A bare + F 2 /2. (7) 

Here, A bare is the bare cosmological constant, which is assumed to be large and negative, and 
F is the magnitude of a four-form field, which can change its value through the nucleation 
of branes. The change of the field strength across the brane is 

AF = ±q, (8) 

where the "charge" q is a constant fixed by the model. 

In order to explain observations, the discrete spectrum of A has to be very dense, with 
separation between adjacent values 

AA < A , (9) 

which in turn requires that the charge q has to be very small. If this is satisfied, analy- 
sis shows that the flat volume distribution (J3J) is quite generic [18]. But once again, the 
exceedingly small charge q required by the model appears to be unnatural. 3 

In an effort to remedy this problem, Bousso and Polchinski (hereafter BP) extended the 
Brown- Teitelboim approach to include multiple four- form fluxes j23|. They considered a 
model in which J different fluxes give rise to a J-dimensional grid of vacua, each labeled by 
a set of integers n a . Each point in the grid corresponds to a vacuum with the flux values 
F a = n a q a and a cosmological constant 



1 J 1 J 

A = A bare + - Fl = A bare + n 2 a ql (10) 



0=1 



This model is particularly interesting because multiple fluxes generally arise in string theory 
compactifications. The model can thus be regarded as a toy model of the string theory 
landscape. BP showed that with J ~ 100, the spectrum of allowed values of A can be 
sufficiently dense, even in the absence of very small parameters, e.g., with |A& are | ~ 1, 
g« ~ 0.1. 



Some ideas on how such small parameters could arise in particle physics have been suggested in 

HEHE1. 
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In the cosmological context, high-energy vacua of the BP grid will drive exponential 
inflationary expansion. The flux configuration in the inflating region can change from one 
point on the grid to the next through bubble nucleation. Bubbles thus nucleate within 
bubbles, and each time this happens the cosmological constant either increases or decreases 
by a discrete amount. This mechanism allows the universe to start off with an arbitrary 
large cosmological constant, and then to diffuse through the BP grid of possible vacua, to 
generate regions with each and every possible cosmological constant, including that which 
we inhabit. 

Our goal in this paper is to study the volume distribution for A in the BP model. In 
particular, we would like to check whether or not this distribution is approximately flat, 
as suggested by the heuristic argument of jjj [Hj. Until recently, such an analysis would 
have been rather problematic, since the calculation of the volume fractions in an eternally 
inflating universe is notoriously ambiguous. The volume of each type of vacuum diverges 
in the limit t — > oo, so in order to calculate probabilities, one has to impose some kind of 
a cutoff. The answer, however, turns out to be rather sensitive to the choice of the cutoff 
procedure. If, for example, the cutoff is imposed on a constant time surface, one gets very 
different distributions depending on one's choice of the time variable t j24[. (For more recent 
discussions, see jiol ufy.) 

Fortunately, a fully gauge- invariant prescription for calculating probabilities has been 
recently introduced in 28]. It has been tried on some simple models and seems to give 
reasonable results. Here we shall apply it to the BP model. 

As BP themselves recognized, their model does not give an accurate quantitative de- 
scription of the string theory landscape. In particular, it does not explain how the sizes 
of compact dimensions get stabilized. This issue was later addressed by Kachru, Kallosh, 
Linde and Trivedi (KKLT) 29], who provided the first example of a metastable string theory 
vacuum with a positive cosmological constant. Apart from the flux contributions in (|10p. the 
vacuum energy in KKLT-type vacua gets contributions from non-perturbative moduli po- 
tentials and from branes. The 4D Newton's constant in these vacua depends on the volume 
of extra dimensions and changes from one vacuum to another. Douglas and collaborators 
FToL l3lL studied the statistics of KKLT-type vacua. Their aim was to find the number of 
vacua with given properties (e.g., with a given value of A) in the landscape. Our goal here 
is more ambitions: we want to find the probability for a given A to be observed. 

In our analysis, we shall disregard all the complications of the KKLT vacua, with the 
hope that the simple BP model captures some of the essential features of the landscape. At 
the end of the paper we shall discuss which aspects of our results can be expected to survive 
in more realistic models. 

We begin in the next section by summarizing the prescription of Ref. J2IJ for calculating 
probabilities. We shall see that the problem reduces to finding the smallest eigenvalue and 
the corresponding eigenvector of a large matrix, whose matrix elements are proportional 
to the transition rates between different vacua. The calculation of the transition rates for 
the BP model is reviewed in Section III. Some of these rates are extremely small, since 
the upward transitions with an increase of A are exponentially suppressed relative to the 
downward transitions. In Section IV we develop a perturbative method for solving our 
eigenvalue problem, using the upward transition rates as small parameters. This method is 
applied to the BP model in Section V. We find that the resulting probability distribution 
has a very irregular, staggered appearance and is very different from the flat distribution 
(JSJ. The implications of our results for the string theory landscape are discussed in Section 
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VI. 



II. PRESCRIPTION FOR PROBABILITIES 



Here we summarize the prescription for calculating the volume distribution proposed 
in Ref. [28]. Suppose we have a theory with a discrete set of vacua, labeled by index j. 



The cosmological constants Aj can be positive, negative, or zero. Transitions between the 
vacua can occur through bubble nucleation. The proposal of Ref. |28[ is that the volume 
distribution is given by 

PjKPjZj, (ii) 

where Pj is the relative abundance of bubbles of type j and Zj is (roughly) the amount of 
slow-roll inflationary expansion inside the bubble after nucleation (so that Zj is the volume 
slow-roll expansion factor). 

The total number of nucleated bubbles of any kind in an eternally inflating universe is 
known to grow without bound, even in a region of finite comoving size. We thus need to 
cut off our count. The proposal of |28[ is that the counting should be done at the future 
boundary of spacetime and should include only bubbles with radii greater than some tiny co- 
moving size e. The limit e — > should then be taken at the end. (An equivalent method for 
calculating pj was suggested in [33].) It was shown in that this prescription is insensitive 
to the choice of the time coordinate and to coordinate transformations at future infinity. 

The bubble abundances pj can be related to the functions fj(t) expressing the fraction 
of comoving volume occupied by vacuum j at time t. These functions obey the evolution 
equation j34| 

|f ij /; • ./ i •/ ' ) • (12) 

i 

where the first term on the right-hand side accounts for loss of comoving volume due to 
bubbles of type i nucleating within those of type j, and the second term reflects the increase 
of comoving volume due to nucleation of type-j bubbles within type-i bubbles. 

The transition rate is defined as the probability per unit time for an observer who is 
currently in vacuum j to find herself in vacuum i. Its magnitude depends on the choice of 
the time variable t. The most convenient choice for our purposes is to use the logarithm of 
the scale factor as the time variable; this is the so-called scale-factor time. With this choice, 

K ij = "^■?~3~^i ' (1^) 

where is the bubble nucleation rate per unit physical spacetime volume (same as \j in 
@) and 

^ = (A,/3)^ (14) 

is the expansion rate in vacuum j. 

We distinguish between the recyclable, non-terminal vacua, with Aj > 0, and the non- 
recyclable, "terminal vacua", for which Aj < 0. Transitions from either a flat spacetime 
(Aj = 0), or a negative A FRW spacetime (Aj < 0), which increase A have a zero probability 
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of occurring. 4 Transitions from Aj = and from small negative Aj to even more negative 
A are possible, but these A's are likely to be large and negative and are therefore of no 
anthropic interest. We will label the recyclable, non-terminal vacua by Greek letters, and 
for non-recyclable, terminal vacua, we will reserve the indices m and n. Then, by definition, 

Latin letters other than m, n will be used to label arbitrary vacua, both recyclable and 
terminal, with the exception of letters a, b, which we use to label the fluxes. 
Eq. (|T2*j) can be written in a vector form, 

f = Mf, (16) 

where f(t) = {fj(t)} and 

Mij — Kij 5ij ^ K r i. (17) 
r 

The asymptotic solution of (jl6j) at large t has the form 

f(t) =f(°) + se- 9 * + ... (18) 

Here, f ^ is a constant vector which has nonzero components only in terminal vacua, 

/f = 0, (19) 

while the values of ffi depend on the choice of initial conditions. It is clear from Eq. ()15|) 
that any such vector is an eigenvector of the matrix M with zero eigenvalue, 

Mf = 0. (20) 

As shown in j2^|, all other eigenvalues of M have a negative real part, so the solution 
approaches a constant at late times. We have denoted by — q the eigenvalue with the 
smallest (by magnitude) negative real part and by s the corresponding eigenvector. 

It has been shown in [28[ that the bubble abundances pj can be expressed in terms of the 
asymptotic solution (fTSj) . The resulting expression is 

(21) 

a 

where the summation is over all recyclable vacua which can directly tunnel to j. 

Note that the calculation of pj requires only knowledge of the components s a for the recy- 
clable vacua. The evolution of the comoving volume fraction in these vacua is independent 
of that in the terminal vacua. Formally, this can be seen from the fact that the transition 
matrix M in (JT5|l has the form 

M -(?o) < 22 > 



4 This is because the volume of the instanton is compact whilst the volume of the Euclideanized background 
spacetime is infinite, so that the difference in their actions is infinite. 
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Here, R is a square matrix with matrix elements between recyclable vacua, while the matrix 
elements of T correspond to transitions from recyclable vacua to terminal ones. Eigenvectors 
of M are of the form f = (s, t), where s is an eigenvector of R, 

Rs = -qs, (23) 

and t is arbitrary. Eigenvalues of M are the same as those of R, except for some additional 
zero eigenvalues with eigenvectors which only have nonzero entries for terminal vacua. 

The problem of calculating pj has thus been reduced to finding the dominant eigenvalue 
q and the corresponding eigenvector s of the recyclable transition matrix R. In the following 
sections we shall apply this prescription to the BP model. 



III. NUCLEATION RATES IN THE BP MODEL 

In the BP model, we have a J-dimensional grid of vacua characterized by the fluxes 
F a = n a q a and vacuum energy densities given by Eq. (JTUJ). BP emphasized that q a need 
not be very small, q a /\-^-bare\ ~ 0.1 — 1. So the model does not have any small parameters, 
except perhaps the values of Aj in some vacua, where the contribution of fluxes is nearly 
balanced by Ab are . 

Transitions between the neighboring vacua, which change one of the integers n a by ±1 can 
occur through bubble nucleation. The bubbles are bounded by thin branes, whose tension 
r a is related to their charge q a as 

rl = ql/2. (24) 

The latter relation is suggested by string theory jijj, 23||. It applies only in the super- 
symmetric limit, but we shall neglect possible corrections due to supersymmetry breaking. 
Transitions with multiple brane nucleation, in which |An a | > 1 or several n a are changed at 
once, are likely to be strongly suppressed [35j, and we shall disregard them here. 
The bubble nucleation rate Ty per unit spacetime volume can be expressed as [36| 

.l„cxp /; (25) 

with 

B id = Iij - Sj (26) 
Here, 1^ is the Coleman-DeLuccia instanton action and 

S i - ~% (27) 

3 

is the background Euclidean action of de Sitter space. 

In the case of a thin-wall bubble, which is appropriate for our problem, the instanton 
action 1^ has been calculated in Refs. (l7l |36| . It depends on the values of A inside and 
outside the bubble and on the brane tension r. 

Let us first consider a bubble which changes the flux a from n a to n a — 1 (n a > 0). The 
resulting change in the cosmological constant is given by 

|AA Q | = K-l/2)g2, (28) 
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and the exponent in the tunneling rate can be expressed as 



'<4 

Here, B^, atspace is the flat space bounce action 



Bai = B flats P a C e r ^ yy (2Q) 



97 2 _4 
yflatspace " ' " 'a 



2 |AA a | 3 

With the aid of Eqs. (j2H) , (j2U) it can be expressed as 

27tt 2 I 



Tjjlatspace 'a ^QH^ 

^al - "IT" | a A I.V ^ dU J 



Tjflatspace " - 1 - /ni \ 

^ 8 K-l/2)3g2 



The gravitational correction factor r(x : y) is given by [3 



rfx v ) = 2 K 1 + ^)-( 1 + 2 ^ + ^] (32) 
x 2 {y 2 -l)(l + 2xy + x 2 ) 1 2 

with the dimensionless parameters 

X = 8jAAJ " 8(n a - 1/2) (33) 

and 

^l^f- 1 - (34) 

where A is the background value prior to nucleation. 
The prefactors Ay in (|2*5j) can be estimated as 

Aij ~ 1. (35) 

This is an obvious guess for nucleation out of vacua with Aj ~ 1. (This guess is supported 
by the detailed analysis in Ref. [Hf.) For Aj <C 1, we still expect Eq. (j33j) to hold, since we 
know that the tunneling rate remains finite in the limit Aj — >• 0, |AA a | ~ 1. 

If the vacuum n a — 1 still has a positive energy density, then an upward transition from 
n a — 1 to n a is also possible. The corresponding transition rate is characterized by the same 
instanton action and the same prefactor j3^|, 

.Zjj Ijii A-ij Aji, (36) 

and it follows from Eqs. ()25j1 . (f27)J) and (fTljl that the upward and downward nucleation rates 
are related by 



Lji — Tjj exp 



24vr 2 



1 1 

a~" a; 



(37) 



The exponential factor on the right-hand side of (|37|) depends very strongly on the value 
of Aj. The closer we are to Aj = 0, the more suppressed are the upward transitions j — > i 
relative to the downward ones. 
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f (A/q 2 ,n) 



200 r 




FIG. 1: The factor f(A/q% function of A/q 2 for n a = 1 (solid line), n a = 2 (dashed line), 

and n a = 10 (dotted line). 

Eq. (|H7|) shows that the transition rate from n a up to n a+ i is suppressed relative to that 
from n a+ i down to n a . It can also be shown that upward transitions from n a to n a+ i are 
similarly suppressed relative to the downward transitions from n a to n a -\. Using Eqs. (J28|) - 
(j23J), the ratio of the corresponding rates can be expressed as 

ln(r i /r T ) = A- 1 /(A/^,n a ). (38) 

The factor f(A/q%, n a ) is plotted in Fig. [I] as a function of Ajq\ for n a — 1 and n a = 2. The 
plot shows that upward transitions are strongly suppressed, unless Ajq\ is very large. The 
factor A" 1 in Eq. (|38|) results in even stronger suppression when A is well below the Planck 
scale. 

Transition rates from a given vacuum j to different states i are related by 

Tij oc exp(— ijj). (39) 

As a rule of thumb, 

h - ~A m L, (40) 

where A max is the larger of Aj and Aj. It follows from (|3*9 ^ .p0 ]l that upward transitions 
from a given site are more probable to the lower-energy vacua. 

To develop some intuition for the dependence of the tunneling exponent B a ^ on the 
parameters of the model, we shall consider the limits of small and large A. For A <C |AA a |, 
we have y ~ — 1, and Eq. (|3~2"j) gives 

r(y -> -1) = (1 -x)- 2 > 1. (41) 

Hence, for low-energy vacua the tunneling exponent is increased over its flat-space value, 
resulting in a suppression of the nucleation rate. (For small values of x, r is increased only 
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FIG. 2: The gravitational factor r as a function of A/|AA | for n a = 1 (solid line), n a = 2 (dashed 
line), and n a = 10 (dotted line). 



The gravitational factor r is plotted in Fig. El as a function of A/|AA | for n a = 1, 
n a = 2 and n a = 10 (corresponding to x = 3/4, 1/4 and 0.04, respectively). We see that 
for large values of A, r <C 1, so the nucleation rate is enhanced. In order for our model to 
be viable, we must ensure that the tunneling action is large enough to justify the use of the 
semi-classical approximation: B a ^ ^> 1, or A/q% <C 20q a ^ 3 . 

If q a and A are changed simultaneously, keeping the ratios A/q% fixed, then x and y do 
not change, and it is clear from Eqs. (|2l?|) . (jHT|) . (jn^)) that the nucleation exponents scale as 
Bij oc A -1 . This shows that bubble nucleation rates are strongly suppressed when the energy 
scales of q a and A are well below the Planck scale. 

IV. PERTURBATION THEORY 

A. Degeneracy factors 

We shall assume for simplicity that the integers n a take values in the range |n | < N, 
where N is independent of a. The number of vacua in the grid is then (2A^ + 1) J . 



by a small fraction, but the factor B a ^ space that it multiplies in the tunneling exponent is 
typically rather large, so the resulting suppression can still be significant.) 
In the opposite limit, A ^> |AA |, y 1, 




and 




(43) 
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To maximize computational abilities, we used the symmetry n a — > —n a and restricted 
the analysis to the sector < {n a } < N. We took into account the degeneracies in A that 
would occur if we allowed negative values of n a by assigning appropriate degeneracy factors 
to the probabilities that we calculated. For example, if we have a two-dimensional grid, 
J = 2, and only consider the quadrant n a > 0, then any point that lies on one of the two 
axes will be doubly degenerate (configuration {0, 1} has the same A as {0, —1}), whilst a 
point that lies in the interior of the quadrant will have a four-fold degeneracy (configuration 
{1, 2} has the same A as {-1, 2}, {-1, -2}, and {1, -2}). 

In general, the degeneracy of each site can be calculated using the following formula: 

V{n a } = 2 k{na} , (44) 

where 

k{n a } = J - (Soni + <W + ••• + S 0nj ) (45) 

So points which have no zero coordinates for a J=7 model have T> = 2 7 = 128. A point with 
one zero coordinate has T> = 2 6 = 64 etc. 

When we use Eq. (J2U) to calculate the probabilities, we multiply the RHS by the appro- 
priate degeneracy. 

Diffusion from a grid point for which n = to n a = —1, is equivalent to the diffusion from 
n = to n a = +1. Also, diffusion from n a = — 1 to n = is equivalent to n a = 1 to n = 0. 
Thus we were able to take into account these processes in the transition matrix by double 
counting the positive n a to or from n = transition rates. In summary, we implemented 
boundary conditions such that our process is equivalent to diffusion through a J-dimensional 
grid, with —N < n a < N. 

As an illustrative example, we show in Fig. El a histogram of the number of vacua vs. A 
for a model with J = 7 and A^ = 4, which has ~ 10 7 vacua. The parameter values used in 
this model are 

q a = 0.5308, 0.3909, 0.5175, 0.4722, 0.5103, 0.4036, 0.4541; A bare = -0.702. (46) 

The sharp spikes and dips in the histogram are due partly to the non-uniform distribution 
of the vacua along the A-axis and partly to the difference in degeneracy factors for different 
vacua. The spikes disappear when the histogram is plotted with a larger bin size, as shown 
in FiglU 



B. Zeroth order 

As outlined in Section II, the calculation of probabilities reduces to finding the smallest 
eigenvalue and the corresponding eigenvector for a huge M X J\f (recycling) transition ma- 
trix R. Here, M is the number of recyclable vacua, which we expect to be comparable to 
the total number of vacua. In our numerical example M ~ 10 7 , while for a realistic string 
theory landscape it can be as large as 10 500 0, B Q IIJ • Matters are further complicated 
by the fact that some of the elements of R are exceedingly small. For example, it follows 
from Eq. (|37|) that upward transitions from low-energy vacua with Aj 1 are very strongly 
suppressed. Matrix diagonalization programs like Mathematica are not well suited for deal- 
ing with such matrices. We shall see, however, that the smallness of the upward transition 
rates can be used to solve our eigenvalue problem via perturbation theory, with the upward 
transition rates playing the role of small expansion parameters. 
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FIG. 3: The spectrum of vacua for a J = 7, N = 4 BP grid with parameters given in (|46[). 



Number of vacua 



800000 - 



600000 



400000 



200000 - 



10 



FIG. 4: The smoothed spectrum for the above model. 



We represent our transition matrix as a sum of an unperturbed matrix and a small 
correction, 

R = R (0) + R (1) , (47) 

where R^ ** contains all the downward transition rates and contains all the upward 
transition rates. We will solve for the zero'th order dominant eigensystem {g( ),s( )} from 
R^ and then find the first order corrections by including contributions from Note 
that the eigenvalue correction is not needed for the calculation of bubble abundances 
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(J21j) to the lowest nonzero order. One only needs to calculate the eigenvector correction s^ 1 ) 

(since the zeroth-order components Sa vanish for recyclable vacua). 
If the vacua are arranged in the order of increasing A, so that 

Ai < A 2 < ... < A^, (48) 

then R(°) is a triangular matrix, with all matrix elements below the diagonal equal to zero. 
Its eigenvalues are simply equal to its diagonal elements, 

flS = -l>i« = -A». (49) 

j<a 

Hence, the magnitude of the smallest zeroth-order eigenvalue is 

g (°) = D am = mm{D a }. (50) 

Up to the coefficient (47r/3)if ~ 4 , D a is the total decay rate of vacuum a. As we discussed 
in Section III, bubble nucleation rates are exponentially suppressed in low-energy vacua with 
Aj <C 1. We therefore expect that the vacuum a* corresponding to the smallest eigenvalue 
g(°) is one of the low-energy vacua. 

With A Qt <C 1 and q a not very small, Eq. (j28|) suggests that downward transitions from 
a* will bring us to the negative-A territory of terminal vacua. Terminal vacua do not belong 
in the matrix R; hence, Rp at = for j3 ^ a*, and it is easy to see that our zeroth order 
eigenvector has a single nonzero component, 

4 0) = <w (si) 

Eq. (j21J) then implies that the only vacua with nonzero probabilities at zeroth order are the 
negative-A descendants which can be reached by a single downward jump from the dominant 
vacuum a*. (Note that the vacuum a* itself has zero probability at this order.) 



C. First order 



The full eigenvalue equation can be written as 

(R(°) + R( 1 ))( S W + s«) = + g( 1 ))(sW + s^). (52) 

Neglecting second-order terms and using the zeroth-order relation 

R (o) s (o) = _ g (o) s (o) j (53) 

we obtain an equation for the first-order corrections, 

(R<°> + g(°)l)s« = -(R« + g«I) S W, (54) 

where I is the unit matrix. 

Eq. (I54j) is a system of J\f linear equations for the J\f components of s'- 1 -'. Note, however, 
that the triangular matrix multiplying on the left-hand side has a zero diagonal element, 

(R(°)W O) I)a,«* = 0, (55) 
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log[l/p 3 ] 




FIG. 5: Plot of log 10 (l/pj) vs. Aj for the BP model with parameters given in (|46l) . The star marks 
the dominant vacuum a*. Triangles represent vacua in group 1, squares in group 2, diamonds in 
groups 3 and 6, crosses in groups 4 and 7, and points in groups 5 and 8. 

which means that the determinant of this matrix vanishes, so it cannot be inverted. In other 
words, the equations in (j54J) are not all linearly independent. 

This problem can be cured by dropping the a*-th equation in (jo^j) and replacing it by a 
constraint equation, which we choose to enforce the orthogonality of and s^ ', 

(sW,s«) = 0. (56) 

Note that the a*-th equation is the only equation in (p>4*j) that involves the eigenvalue cor- 
rection Now has dropped out of our modified system, and we can straightforwardly 
solve for s« . We did this numerically for a J = 7 model; the results will be presented in 
the following subsection. A J = 2 analytic toy model is worked out in the Appendix. 



V. BUBBLE ABUNDANCES IN THE BP MODEL 

We found in the preceding section that the zeroth order of perturbation theory picks the 
vacuum a* which decays the slowest (we call it the dominant vacuum), and assigns non-zero 
probabilities to it's offspring only - all other probabilities are zero. In the first order of 
perturbation theory, all vacua connected to the dominant vacuum via one upward jump, 
and any vacua connected to these via a series of downward transitions, also acquire non-zero 
probabilities. 

The bubble abundance factors pj for the 7-dimensional toy model (J4~rj|) are shown in Fig. 
El We plot log 10 (l/pj) vs. Aj, so higher points in the figure correspond to smaller bubble 
abundances. The first thing one notices is that there are several groups of points, marked 
by triangles, boxes, etc. The star marks the dominant vacuum a*. 
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In this particular example, the dominant site has coordinates (1,1,1,1,1,1,1). There 
are J = 7 ways to jump up one unit from this site, to arrive at the seven points indi- 
cated by black triangles, which we shall call group 1. The coordinates of these points are 
(2,1,1,1,1,1,1), (1,2,1,1,1,1,1), (1,1,1,1,1,1,2). Lower-energy vacua in this group 
have higher bubble abundances, in accordance with Eqs. (J3T?|) . (J4TH) of Section III. 

The next group of states results from downward jumps out of vacua in group 1 in all 
possible directions, excluding the jumps back to the dominant site a*. We call it group 2. 
The number of states in this group is J(J — 1) = 42. Consider, for example, the subgroup 
of states in group 2 coming from the downward transitions out of the state (2, 1, 1, 1, 1, 1, 1). 
These states have coordinates (2, 0, 1, 1, 1, 1, 1), (2, 1, 1, 1, 1, 1, 0). Since they originate 
from the same single parent, the difference in their bubble abundances comes from the 
difference in the instanton actions 7y. This effect is much milder for downward transitions 
than it is for the upward ones. That is why the spread in bubble abundances within the 
subgroups of group 2 is much smaller than it is in group 1. 

Further downward jumps replacing one of the J — 2 = 5 remaining l's by a give rise to 
group 3, consisting of J(J — 1)(J — 2)/2! = 105 states having flux configurations with one 
n = 2, fourn = 1 and two n = 0. Similarly, group 4 includes J(J — I) (J — 2)(J — 3)/3! = 140 
states with one n = 2, three n = 1 and three n = 0, and group 5 includes J (J — 1)(J — 
2) (J — 3) (J — 4)/4! = 105 states with one n = 2, two n = 1 and four n = 0. The factorial 
factors are included to avoid double counting. For example, the site (2, 0, 1, 1, 0, 1, 1) can be 
reached by downward jumps from either (2, 1, 1, 1, 0, 1, 1) or (2, 0, 1, 1, 1, 1, 1) and would be 
counted twice if we did not divide by 2!. 

If a vacuum in group 2 has a coordinate jump from n = 2 to n = 1, the resulting site is 
one of the daughter sites which can also be reached by downward jumps from the dominant 
site. These negative-A vacua have non-zero probabilities already at the zeroth-order level 
and are not represented in the figure. 

If a vacuum in group 3, 4 or 5 has a coordinate jump from n — 2 to n — 1, the resulting 
sites are all terminal vacua (groups 6, 7 and 8, respectively). 

We note that although the dominant vacuum a* is one of the low-energy vacua, there 
are many other recyclable vacua which have lower A. Recall that the dominant vacuum 
has the smallest, in magnitude, sum of its transition rates down in each possible direction 
(see Eq.'s fl4T?|) and (jSOI))- Each transition rate depends exponentially on the value of q a 
and the factor r(x,y)/(n a — 1/2) 3 . From this factor and Fig.Q we see that for A/AA < 1 
(this is the case for a*) any jump from an n = 2 flux quanta will be less suppressed than 
a jump in the same direction from an n = 1 flux quanta. Thus we are not surprised that 
states which contain a flux quanta of n = 2 are not dominant sites despite having smaller A 
than a*. Since each transition rate is exponentially dependent on the tunneling exponent, 
typically the largest (in magnitude) transition rate will dominate the sum in Eq.()49|). Thus, 
essentially for a vacuum to be the dominant state its largest (in magnitude) transition rate 
should be smaller that the largest transition rate of any other vacuum. 

The distribution in Fig. |5]was obtained in the first order of perturbation theory, which 
includes only the vacua which can be reached by a single upward jump from the dominant 
site a*, followed by some downward jumps. If higher orders were included, we would see 
additional groups of vacua, reachable only with two or more upward jumps. These vacua 
would have much smaller bubble abundances than those already in the figure. 

The distribution pj in Fig. El spans more than 300 orders of magnitude. It differs dramat- 
ically from the flat distribution (jSJ) suggested by the heuristic argument in the Introduction. 
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log[l/ Pj ] 




FIG. 6: Plot of log 10 (l/pj) vs. Aj for the BP model with parameters (|57jl. The star marks the 
dominant vacuum a*. Different groups of vacua are represented by the same symbols as in Fig. 5. 

Many vacua with close values of Aj have very different abundances pj. The reason is that 
despite their closeness in A, such vacua are located far away from one another in the BP 
grid, and the paths leading to them from the dominant vacuum a* are characterized by 
exponentially different transition rates. Even the vacua resulting from tunneling out of the 
same site typically have very different abundances, due to the exponential dependence of 
the tunneling rates on q a . 

Fig. IH1 shows the distribution of bubble abundances for another J = 7 model, with a 
different set of parameters: 

ft = 0.6175, 0.3909, 0.6472, 0.5508, 0.5103, 0.7036,0.4541; A bore = -1.033. (57) 

In this case, there is more scatter in the values of q a , and the groups of vacua are somewhat 
less pronounced. However, the staggered nature of the distribution is still apparent. 



VI. DISCUSSION 

In this paper we have used the prescription of Refs. jH, |33| to determine the bubble 
abundances pj in the BP model. We found that the resulting distribution is very irregular, 
with values of pj soaring and plummeting wildly as Aj changes from one value to the next. 
This distribution is drastically different from the flat distribution (J5J which was used as a 
basis for the successful anthropic prediction for A. 

Apart from the bubble abundance factor pj, the volume distribution (fTTj) includes the 
slow-roll expansion factor Zj. In any realistic model, bubble nucleation should be followed by 
a period of slow-roll inflation, at least in some bubble types. The expansion factor Zj is, of 
course, model-dependent, but there is no reason to expect that it will somehow compensate 
for the wild swings in the values of pj as we go from one value of Aj to the next. 
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Another point to keep in mind is that, in a realistic setting, vacua with different values 
of the fluxes F a may have different low-energy physics, so the density of observers n & s (A) 
would also be very different. We should therefore focus on the subset of vacua in the BP grid 
which differ only by the value of A and have essentially identical low-energy constants. Once 
again, there seems to be no reason to expect any correlation between these constants and the 
up and down swings in the bubble abundances. We conclude that the staggered character 
of the distribution Pj = P(Aj) is expected to persist, even in more realistic versions of the 
model. 

This conclusion is not limited to the BP model. It is likely to arise in any landscape sce- 
nario, where a dense spectrum of low-energy constants is generated from a wide distribution 
of states in the parameter space of the fundamental theory. Vacua with nearly identical 
values of A may then come from widely separated parts of the landscape and may have very 
different bubble abundances and volume fractions. 

Given the staggered character of the volume distribution, what kind of prediction can 
we expect for the observed value of A? The answer depends on the number Ma of possible 
vacua with Aj within the anthropic range (jljl. AA^ ~ 10~ 118 . (We count only vacua in 
which all low-energy constants other than A have nearly the same values as in our vacuum.) 

Suppose the volume factors in the distribution Pj span K orders of magnitude. (K ~ 300 
in our numerical example in Section V.) We can divide all vacua into, say, 10K bins, such 
that the values of Pj in each bin differ by no more than 10%. Suppose now that there are 
Ma "C 10 K vacua in the anthropic range AA^. We can then expect that most of these vacua 
will be characterized by vastly different volume factors Pj, so that the entire range will be 
dominated by one or few values of Aj having much higher volume fractions than the rest. 

Moreover, there is a high likelihood of finding still greater volume fractions if we go 
somewhat beyond the anthropic range - simply because we would then search in a wider 
interval of A. We could, for example, find that a vacuum with Ai ~ 10~ 114 ~ 10 6 Ao 
has a volume fraction 200 orders of magnitude greater than all other vacua in the range 
< A < Ai. Galaxy formation is strongly suppressed in this vacuum: the fraction of matter 
that ends up in galaxies is only fa(Ai) ~ 10~ 110 . However, this suppression is more than 
compensated for by the enhancement in the volume fraction. 

If this were the typical situation, most observers would find themselves in rare, isolated 
galaxies, surrounded by nearly empty space, all the way to the horizon. This is clearly not 
what we observe. The dominant value could by chance be very close to A = 0, but if such 
an "accident" is required to explain the data, the anthropic model loses much of its appeal. 

In the opposite limit, Ma ^> 10K, the number of vacua in the anthropic interval AA^ 
is so large that they may scan the entire range of Pj many times. Then, it is conceivable 
that the distribution will become smooth after averaging over some suitable scale SA. If 
5A can be chosen much smaller than AA^, then it is possible that the effective, averaged 
distribution will be flat, as suggested by the heuristic argument in the Introduction. The 
successful prediction for A would then be unaffected. 5 

The above argument is somewhat simplistic, as it assumes that the vacua in the BP grid 
are more or less randomly distributed between the 10K bins, with roughly the same number 
of vacua in each bin. Such an "equipartition" is not likely to apply to the most abundant 



5 Joe Polchinski has informed us that a similar argument, indicating that the anthropic explanation for 
the observed A requires a large number of vacua in the anthropic range, was suggested to him by Paul 
Stcinhardt. 
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vacua, but it may hold for the vacua in the mid-range of Pj. Finding the conditions under 
which equipartition applies would require a statistical analysis that goes beyond the scope 
of the present paper. 

In summary, it appears that the staggered volume distribution resulting from the BP 
model is in conflict with observations, unless it yields a huge number of vacua in the anthropic 
range of A. Counting only vacua which have nearly the same low-energy physics as ours, we 
should have much more than 10K; hence, the total number of vacua should be many orders 
of magnitude greater. The large number of vacua in the anthropic range is only a necessary 
condition for the distribution Pj to average out to the flat distribution (jSJ). Further analysis 
will be needed to find whether or not this actually happens, and if so, then under what 
conditions. It would also be interesting to analyze other simple models of the landscape, 
such as the "predictive landscape" of Arkani-Hamed, Dimopoulos and Kachru 44l |. and see 
what similarities and differences they have compared to the BP model. 

Throughout this paper we assumed that the brane charges q a are not particularly small. 
This assumption may be violated in certain parts of the landscape, e.g., in the vicinity 
of conifold points, resulting in a much denser spectrum of vacua I20L l3ll . I32L |41 |. Infinite 
accumulations of vacua may occur near certain attractor points [42|, |43[. Implications of 
these effects for the probability distribution on the landscape remain to be explored. 
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VIII. APPENDIX 

We will illustrate our perturbative method of calculation on a very basic BP model, 
which can be solved analytically. We consider 9 vacua arranged in a 2-D grid and labeled 
as indicated in Fig. [7| There are three terminal vacua labeled 1, 2, 4, and six non-terminal 
vacua, 3, 5, 6, 7, 8, 9 in this model. We allow upward and downward transitions between 
nearest neighbor pairs, with transitions from non-terminal to terminal states allowed, but 
no transitions may take place from a terminal state. For simplicity, we disregard the vacua 
in the quadrants where rii < and/or n 2 < and assume that the set of vacua in Fig.(J7J) is 
all there is. 

The evolution equations for the set of non-terminal vacua is 

| = Rf, (58) 

where the vector f = {/ 3 , / 5 , / 6 , fr, f$, fg}. Assuming that upward transition rates are far 
more suppressed than downward transition rates, we represent the transition matrix as 

R = R(°) + RW (59) 
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FIG. 7: The arrangement of vacua for a J = 2, N = 2 BP grid 

where 



R(°> = 

and 
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and we have defined 

D a = J2 Kja, (62) 

j<a 
j>a 

D a and U a represent, respectively, the sum of the downward and upward transition rates 
from vacuum a. 
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In our toy model we will assume that vacuum 5 has the smallest (in magnitude) sum 



of downward transition rates, and therefore — q 



(0) 



-Dk is the zero'th order dominant 



eigenvalue. By inspection, we see that the corresponding eigenvector is = {0, 1, 0, 0, 0, 0}. 

We now need to calculate the first order correction to this eigenvector, s^. Substituting 
(|dTH1 . (|dT|) in Eq. we find 
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Note that the only equation in this set that depends on the first order correction to the 
eigenvalue is also the equation that needs to be dropped from our system, since it has a 
zero diagonal element - this causes the matrix on the right-hand side of (JMj) to have a zero 
determinant, which renders it uninvertible. 

This drop in the number of independent equations is replenished by including the con- 
straint equation (jo^j) ; the resulting set of equations is 
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The solution is readily determined, and we obtain 
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This can now be used in Eq. fl21|) to determine the bubble abundances. For example, com- 
paring the bubble abundances in vacua 3 and 7, we find 



P3 _ HqK 36 k 65 (D s — D 5 ) 
p 7 HlK 1& K Sb (D & - D 5 ) 



(67) 
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